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SUMMARY 


. A methodology for simulating the growth of long through cracks in the skin of pressurized 

^pSt o f f ?SSr e Th described. Crack trajectories are allowed to be arbitrary and are computed 
as part of the simulation. The interaction between the mechanical loads acting on the superstructure and 

strategy^ S SSctoTre S S n near ? e Cra n k ti P s , is ' accounted for by employing a hierarchical modeling 
strategy. 1 he structural response for each cracked configuration is obtained using a geometricallv & 

non inear shell finite element analysis procedure. Four stress intensity factors? I ^f^membrnne 
behavior and two for bending using Kirchhoff plate theory, are computed using an extension of the 
modified crack closure integral method. Crack trajectories are deteSed by fptf yh^ 

renSSied Sti f eSS c 5 lte y ion - Crack growth results in localized mesh deletion, and the deletion regions are 
ofThp^ ^ih U ?? atlCa11 ^ USmg ^ newly developed all-quadrilateral meshing algorithm. The effectiveness 
SemonSed^ appllC ? bllit y to Performing practical analyses of realistic structuresls 
ty w^ZZ ° c ^ a ek growth in a fuselage panel that is representative of a 

yp cal narrow-body aircraft. The predicted crack trajectory and fatieue life comnare well with 
measurements of these same quantities from a full-scale pressurized paneltest 


INTRODUCTION 


i Coocern for the safety of aging airliners motivates the development of a methodolnpv fnr 
simulaung the behavior of long through cracks in the skin of piessuSffi 

Zrinfffie^ Long cracks in the fuselage skin can arise from accidental damage 

ng flight, e.g., the failure of an engine component resulting in a puncture of the skin or from fatimip 
crack growth whereby several smaller cracks along a row of n* ets SVtove linked^’ TlSfn n? a 

ForS e ri a m Ia m S p ? ssu ? 2ed - Since the loadi "S from pressurization is cyclic 1 ^^ 0 ^^^^^^^;^^ 

grow* Re?mdtess r of fte cause ScTokf 5 sig , nifiCant n T- bers ° f Cycles and substantial fatigue crack 
aircraft d^ag^tolerarice!^ 0 ^ 0 ^ “““ he ' P qUantify the effect ot ma j° r desi «" Parameter changes on 
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Crack growth simulations, in general, are difficult to perform using existing codes for 
computational solid mechanics alone. The difficulty lies in the need for a continual updating of the 
geometry, as well as the discretization of the structure, as a result of crack initiation and propagation. 
Most available codes will not allow the geometry of the structure to evolve automatically during the 
course of a simulation, thereby forcing the analyst to perform tedious and time-consuming manual 
updates to the model for each step of crack growth. These difficulties can be overcome by use of a 
single fracture simulation program containing a representational model of the problem that separates the 
problem description from the problem discretization. Such a program provides the necessary framework 
for the entire simulation process, employing an existing computational solid mechanics module simply 
as an analysis engine whenever the structural response of a particular cracked configuration is necessary. 
The representational model consists of a geometrical model to which simulation attributes, such as loads 
and boundary conditions, can be attached once at the start of the simulation. In this way , information is 
inherited automatically by each computational mesh that is generated during crack growth. The primary 
components of this discrete crack growth analysis methodology are shown in Figure 1. It has been 
implemented within the FRANC3D fracture simulation program [Wawrzynek, et al. (1988); Martha 
(1989)- Wawrzynek (1991)], a general software framework that supports simulation of curvilinear 
discrete crack propagation in stiffened, thin shell structures [Potyondy (1993)] and non-planar discrete 
crack propagation in solid structures [Potyondy, et al. (1992); Martha, et al. (1993); Gray, et al. (1994)]. 


The discrete crack growth analysis methodology described here utilizes the FRANC3D software 
framework. It accounts for the dominant influence factors affecting crack growth behavior and 
integrates these with a set of fracture parameters that characterize the fracture process. The mechanics 
considerations underlying the methodology are summarized in the next section. This is followed by a 
description of the modeling approach that supports simulating crack growth in realistic fuselage panels. 
Then, two example simulations are presented in which computed stress intensity factor histories are used 
to predict crack trajectory and fatigue life and a comparison is made with measurements from ^ 
scale pressurized panel test. A more thorough description of this work can be found in Potyondy (1993) 
and Potyondy, et al. (1994b). 


MECHANICS CONSIDERATIONS 


The behavior of long through cracks in the skin of pressurized fuselages is dominated by the 
following four factors [Rankin, et al. (1993)]: 1) & geometrically nonlinear stiffening effect (hereafter 
referred to as the bulge-out effect), caused by tensile membrane stresses in the direction parallel to the 
crack that restricts the crack edge bulging; 2) the mechanical loads, in addition to the internal pressure, 
acting on the superstructure that alter the stress distribution in the skin; 3) the presence of stiffening 
elements (frames and stringers) that alter locally the stress distribution in the skin; and, 4) plasticity 
effects around the crack tip. The first three factors can be modeled using a geometrically nonlinear shell 
finite element analysis procedure in which cracks are represented discretely in the mesh. In the present 
work, the nonlinear shell finite element code STAGS [Almroth, et al. (1986)], is used along with a 
hierarchical modeling strategy to account for the second factor. The effect of plasticity at the crack tips 
is ignored. (This effect is accounted for in the work of Shenoy, et al. (1994), who compute T*, a 
relevant no nlin ear fracture parameter, for typical fuselage structures. Refer to that paper for further 
details.) 


The use of a geometrically nonlinear finite element shell model captures the global structural 
response by idealizing the fuselage as a cracked, stiffened, thin shell structure; however, additional 
modeling idealizations are required to predict crack trajectory and fatigue life. For cases in which the 
linear thin shell equations are applicable, stress intensity factors can be defined [Potyondy, et al. 



in *2 ' i H°^ ev . er ’ in the present case, wherein geometrically nonlinear effects plav a significant role 

m the shell behavior, there is no clear consensus about the proper definition of s 

the solu^on Droc^nre ln pir^ fiel ? ha ^ e . adopted an engineering approach that restricts the generality of 
. P .st, only straight cracks, oriented and loaded symmetrically with respect to the 

t ^ T Sldere< ? ; thUS ’ a Sin « le similitude parameter, the mode-I stress intensity f altar K can 
crack^JS C t 0rrel . ate crack growth rate with experimentally obtained fatigue crack growth data Also 
effect k C ° nSldered ’ sm . ce die crack growth is assumed to be self-sinhlar. The bulge-out 

effect modeled by using a geometrically nonlinear solution procedure, and the energy release fate G, is 

2^' Kl iS obtained f rom the linear elastic, plane stress relation K r = ^EG where E is Young’s 


k oontrasf the present methodology addresses the more general problem in which the cracks 
may be curved and need not be oriented or loaded symmetrically. Such cracks experience a combination 

and out-of-plane benLg loading" The 

* d case is snown m figure 8* The unsymmetnc out-of-plane bending could be called a nominal 
mode-ni deformatton. Hui and Zehnder (1993) propose a set of foSs Sn^ty teofs? ^ 

,^2 , Kj , Kjj ) to account for this general loading. The mechanical behavior of the region near the 

plane s'tfes^ plate * eoi y J° accoaat for *e bending loads and two-dimensional 

nlatc thco^ H ^ u T 5 , the membrane loa ds- The stress intensity factors based on Kirchhoff 
p te theory are denoted by k x and while those based on two-dimensional elasticity theory are 

dea0 ^ db f K J. and K u- A11 / our stress intensity factors have the units of force/length^ . The four 
is characterized g^e cteptcted in Figure 2. Note that the unsymmetric out-of-plane bending in Figure 8 

^ ^ be stress intensity factors described above arise from the assumptions of small deflection nlate 

geo " ly r nlinear she11 finite elemenf p roced ^ Sffi 

J? } effect experienced by a long crack in a pressurized fuselage structure. In the present work the 
engineering assumption is made that these stress intensity factors can also be used to characterize the 

S&trr i ThUS ’ ““ ***** fact0rs co ^ e ^(m^o°l^^dlibrium 1 state 

the . geometrically nonlinear shell computation. They are computed usina an extension of 
the modified crack closure integral method [Potyondy (1993); Potyondy et al (1993)1 that is 
implemented as a post-processing step to the she'll finite element a"na^ 

element 16185 ^ 6 6aS< : rate *?* f* 6 obtauied from the nodal forces and displacements in the shell finite 
the expSonsf" “ P ' TheSe com P onents «"• in turn, related to the stress intensity factors via 




7= 1 


1 E 

G 2 =^ 

2 E 

n k x jvf 1+ v 

G4= ^“l3^J 


(1) 
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Cr3+Cr6- 3 E U+v 


The stress intensity factors so obtained correspond to an equivalent model of a cracked Kirchhoff plate 
with the same energy release rate as that computed for the geometrically nonlinear problem. 


No generally accepted criterion exists for predicting trajectories of cracks in pressurized thin 
shells. Existing criteria strictly apply for pure membrane stress conditions only. The stress state about a 
crack in a pressurized thin shell will always contain out-of-plane bending components, and additionally, 
these bending components will become large as the crack approaches stiffening elements. It is no ye 
apparent in what way these bending components are to be incorporated into an appropriate criterion. 
Rankin, et al. (1993) comment on this problem and refer to the summary of Zaal (1992), however, no 
suggestions are made for incorporating the bending components. 


In this work, the crack front fields are described by , K u ) . These fields apply to a 

crack in a plate wherein the crack front is assumed to be straight and aligned with the plate normal 
direction. (If the crack front were to deviate from this shape and orientation, a three-dimensional model 

of the crack front stress field involving (Kj ( s ), K u ( s),K m (s)) would be required where s 
parameterizes the crack front through the thickness.) Therefore, crack trajectory can be described u y 
by the path of a single point on the crack front at the shell midsurface. The path is described in terms of 
the instantaneous angle of crack propagation (j), with which the crack will deviate from its current se 
similar direction. $ will lie in the tangent plane to the midsurface at the crack tip. The maximum 
tangential stress theory [Erdogan and Sih (1963)] is used to predict <J) via the relation 


Kj sin (p + K n (3 cos <p - 1) = 0 


( 2 ) 


Since cj) is determined only at the shell midsurface, the two bending stress intensity factors do not enter 
into the relation. 


MODELING APPROACH 


Hierarchical Modeling Strategy 


A hierarchical modeling strategy [Starnes and Britt (1991)] consisting of a series of models at 
differing scales provides a general mechanism for modeling pressurized fuselage structures by 
representing the important global and local features of a shell structure subjected to combined interna 
pressure and mechanical loads. The hierarchical modeling levels range from a relatively ^ coarse global 
shell model for determining the internal load distribution and global response of the shell, to a hig y 
refined local shell model for determining crack growth behavior. 

The three hierarchical modeling levels employed in both of the simulations are depicted in Figure 
3 and include a global shell model, a 6x6 bay stiffened panel model, and a 2x2 bay stiffened panel 
model. The global shell model can be subjected to internal pressure and mechanical loads that represent 
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a bending moment, a vertical shear load, and a torsional load applied to the shell A seometricallv 
nonlinear analysis of the global shell model provides the internal load distribution f of the shell as well as 

^^^el b i 1 IJw y T C01 !M? ?nS (dl t plaC ? ments and rotati ons) along the boundaries of the 6x6 bay 1 
SSt addltl o n to these kinematic constraints, the 6x6 bay stiffened panel model is 

mnl S | b eC ^ t0 mternal Pressure. A geometrically nonlinear analysis of the 6x6 bay stiffened panel 
n P I° v i? es a n ^ or f accura t e determination of the interaction between skin, frames and stiffeners as 

addkfon m SeST’ 110 b °? ndafy oonditions along the boundaries of the 2x2 bay stiffened panel model. In 
tion to these kinematic constraints, the 2x2 bay stiffened panel model is also subiected to internal 
pressure. A geometrically nonlinear analysis of the 2x2 bay stiffened panel model provides the detailed 

lenPth? d A def0r ? a , ti0n ? eld near the crack tips used t0 compute (h>h,Z Zj>K n ) for successive crack 
S S 9v9 ut n C I aCk ^ pres , e ? t in a11 three of the hierarchical models, but the crack is grown only in 

crack ^owth ^tinf the" T ; US ’ il i kin ® matic bounda ry conditions are not recomputed during 
crack growth. All of the analyses are conducted using the STAGS shell finite element analysis code 


Modeling Crack Growth 


FR AN(An C anH k ch r ° Wtk t S 6S pla , ce i? tiie 2x2 bay stiffened panel model that is constructed within 
IC\NC3p and shown m Figure 4. The general steps of the process are illustrated in Figure 1 (In this 

panel t ^d 1 <^an^are atoche^ffTth 11 ^ ^ disp ! acemen ^ obtained from the 6x6 bay stffened 

penmptrhS 1 adached t0 the corresponding edges along the boundary of the FRANC3D 

SZsS'eS the m,alySt Can “ nipulate a view 


process Znsists of itolZ‘fhro n „ a ifi b T. S , haye bee " atta ^ 1 ? 1 , 10 the S eom ettical model, the simulation 
t ? ugh the following steps: 1) A STAGS finite element model is generated 

from FRANC3D for the current cracked configuration. 2) After obtaining the equilibrium state by 

growth diction ^Th^ 6 ’ t 1 he f ValuCS 0f (^>^2 > K i,K u ) are computed and are used to estimate crack 
^ i i ^ The analyst determines the amount of crack extension at each crack tin 4^ The 

SlefioTTlVh 1 ^^, 6 tOP ° 10gy ? f * he ^C3D geometrical model Std lead“m loSnSsf 6 
deletion. 5) The deletion region is remeshed automatically using a newly developed all-auadrilateral 

^r^ CShmg ^ g ° tithm [ p0t y°« d y> «al. (1994a)] which geLatL TgradS 

element fRanldn and^mpa^^w! 1 !! ^ crack .ti ps - A four-noded, six degree-of-freedom per node, shell 
element LKanJan and Brogan (1991)] is used within the entire deletion region including the area 

immedrateiy ad^cent to the crack tips. (The errors introduced in the strefs intensity faftor^ computation 
from elemen distortion are examined in Potyondy (1993). For the meshes us^ Ke ex^pl7 

may reacri5 neSnn 1 Tif 7 ^ K " beUev u ed to be below 2 percent while the errors in k x and \ 

andVZ/s *• Proceeding in this fasUon ' both ** ***** 


Generic Narrow-Body Panel 


.• A tiiselage panel that is representative of a typical narrow-body aircraft has been modeled The 

Ses°s n 0 S 0?6 tST° n r fr ,T MUler ’ * aL (1992 >- The panel consi^f ^20 2 Tm?ad sMn Of 

Spacing tZ sSo^Sn&n"#^ 5 V' 25 inCh and frames - shear-fies, at 20 

P g. 1 ear straps of thickness 0.036 inches run circumferentially at 10 inch spacing, centered on 
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the frame and midway between frames, and longitudinally under the stringers. 


The panel is comprised of the following structural components: skin, frames, stringers, tear 
straps, and stringer clips. These are shown in Figure 4. Typically, the tear straps are bonded dkectly to 
the skin whereas the other components are connected with rivets. In the work of Miller, et al. (1992) 
fastener flexibility was modeled using rigid body elements for rotational and axial displacements, and 
spring elements for transverse displacements. In the present work, fastener flexibility has not been 
modeled; it is assumed that there is a perfect bond between all structural components along then entire 
area of overlap. 


Different modeling idealizations have been employed in the three models of the analysis 
hierarchy. In the global shell model, the skin is modeled with shell elements; the frames are modeled 
with shell elements to which beam elements are attached along some of the shell element edges; and, the 
stringers, tear straps, and stringer clips are modeled with beam elements. In the 6x6 bay stiffened panel 
model, all structural components are modeled with shell elements and beam elements are also included 
along the tops of the stringers. In the 2x2 bay stiffened panel model, these beam elements are replaced 
by shell elements that approximate the hooked shape. In both the 6x6 and 2x2 bay stiffened panel 
models, each layer in the overlap areas is modeled with its own separate shell elements by providing the 
eccentricity from a reference wall surface. The material properties used in the three models of the 
analysis hierarchy were a Poisson’s ratio of 0,33 for all structural components; and, a Young s modulus 
of 10,500 ksi for the skin and tear straps, and 10,700 ksi for the frames, stringers, and stringer clips. 


EXAMPLE SIMULATIONS 


Simulation 1: Trial Problem 


The initial crack is located in a panel near the crown of the fuselage. It is. a straight 6.0 inch, 
crack oriented parallel to the stringer and centered on the midbay tear strap 2.31 inches from the stringer 
centerline. The crack cuts completely through both the tear strap and the skin. The boundary conditions 
for the global shell model are a uniform internal pressure of 8.0 psi and a mechanical load consisting of a 
bending moment applied to the fuselage barrel. A plane of symmetry is assumed at one end of the 
barrel. The other end is restrained to remain in a plane and two equal but opposite point forces are 

applied at the top and bottom of the barrel. The point forces produce a moment of 5. 65 X 10 
pounds, directed so as to compress the fuselage crown. The applied moment is slightly less man the 
value that just causes buckling of the skin between frames in the fuselage crown. The crack is grown in 
one inch increments at each crack tip up to a total piecewise crack length of 16 inches. Both the gL*ba 
shell and the 6x6 bay stiffened panel models contain approximately 100,000 degrees of freedom. 1 he 
FRANC3D models contain from 11,000 to 14,000 degrees of freedom. 


The deformed structure for various crack lengths is shown in Figure 5. The computed crack 
trajectory, shown in Figure 6, appears reasonable. The crack is curving away from the nearest stringer. 
This is the observed behavior in pressurized panel tests of narrow-body fuselage structures wherein the 
crack is located closer to the stringer. (See the description of the test results in a subsequent section and 
Figure 9). 


586 



crack length!?!® tip is Plotted in Figure 7 in terms of relative 

frames. The stress intensity factor historv fnr the. i f t we ® n crack Ups and B is the distance between 
stress intensity factor histo^ can be identified.^ ^ UP 1S SUTular - 1116 following trends in the 

Kj is the dominant stress intensity factor and increases as the crack grows. 

' toou^ou^^ s ^r C mr^ ter * h * f* Ck “ aU ° Wed to CUIVe Md small 

subjected 'to^mpl^el^din^i^^WawrajmeTifl 99 1)] 001 ^ fOT “ 

ra , ti0 Vf' 1S smaU for 1116 initial straight crack, but as the crack begins to curve 
the value of k, mcreases such that the ratio of k^/K, becomes large. 

The ratio k^/Kj remains small throughout crack growth. 

deform^hm^Twobends ^e^xuectokto ° f °PP° sin S. <™ck edges to out-of-plane 

the stiffness of the crack edge nearest to the strip ?e^h a o l0 iH gltUdinal cr f cl f is moved closer to a stringer, 
crack edge; and 2) as the era Hr hpenn* xn crease relative to that of the opposing 

Should dfmirtish relative lo" that of die cracked p’e^n Tr™ ° f *5 0n *» co ™ave ^ 

the longitudinal direction that restrict the buleinp ranJjnt,? 1 ™ 6 !* S1 . de ‘ ! ensile membrane stiesses in 

concave side.) The initial crack in the trial develo P m th L e peninsula of material on the 

to be observed; however, it is observed in the valhStinnn^f en ° ugh J° the stringer for the first trend 
observed in both the deformed shapes of Figure 5 anri SC6 Flg ? re The second tiend is 

as the crack curves, both th^ unsv^eMc J T h ^ faCt ° r hiSt ° ry of Figure 7 > i e ’ 

low valle nfk 1 Th 6 Symn ; etric bending occurring for the crack in the trial problem as indicated bv the 


Simulation 2: Validation Problem 


Description of Experiment 


Group. The fixtures ?an be u£d for E™' fabricated by the Boeing Commercial Airplane 
has a radius of curvature of 74 tacte to mSch ” resid “ al e Strength tests - 0" e fixture 

curvature of 127 inches to nmch wMe“bSv^ J?* the 0ther has a radius of 

framework and the other on rollers to f!r G ^ nd of each flxture 1S mounted in a rigid 

inserted in cutouts in the barrel-shaDed fixture^ T U< ^ na * dlSp a . Ceme , nt ‘ Rem ovable test sections are 
air as the pressurizing medium More informant 6 k‘ S a f®. CO o ducted under pressure loading only, using 
found in the papers o! Miller, ut the )fix ture s and the Boeing test program caTbe 
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Experimental data generated using the test fixtures was made available to the authors by Boeing 
[Worden (1993)]. The data is from a typical narrow-body test panel with bonded tear straps. The panel 
had initial damage consisting of a 5.0 inch saw cut centered on the midbay tear strap and just above the 
stringer tear strap. The saw cut went completely through both the skin and the midbay tear strap. The 
panel was inserted into the test fixture and pressure cycled at 7.8 psi. The pressure cycles fluctuated 
between maximum pressure and zero pressure. During the test, the position of each crack tip was 
monitored and recorded along with the number of pressure cycles. The experimental data can be found 
in Potyondy (1993) and Potyondy, etal. (1994b). 


Description of Simulation 


The test fixture approximates the behavior of an actual pressurized fuselage. Thus, a cracked 
pressurized fuselage is modeled using the analysis methodology described above. The validation 
problem differs from the trial problem only in the following details: 1) initial crack size and location; 
and, 2) boundary conditions applied to the global shell model. The initial straight 5.0 inch crack is 
located in a panel that is representative of the fuselage crown. In both the global shell model and the 
6x6 bay stiffened panel model, the initial crack is located in the same location as the saw cut. In the 
FRANC3D model, the initial crack has been translated a distance of 0.45 inches (half the stringer width) 
away from the stringer, because the current version of FRANC3D cannot model cracks at the interface 
of the stringer and skin. 


The boundary conditions applied to the global shell model are a uniform internal pressure of 7.8 
psi and no mechanical loading. A plane of symmetry is assumed at one end of the barrel while the other 
end is constrained to remain in a plane. Note that this model corresponds with an open cylinder, and 
thus the longitudinal stress in the model is less than the longitudinal stress in the test fixture, since the 
test fixture is a closed cylinder. The crack is grown in one inch increments at each crack tip up to a total 
piecewise crack length of 19 inches. 


The deformed structure for various crack lengths is shown in Figure 8. The computed crack 
trajectory is plotted in Figure 9 along with the measured trajectory from the panel test. The computed 
stress intensity factor history for the right crack tip is plotted in Figure 10. The stress intensity factor 
history at the left crack tip is similar. 


There is no analytical solution with which the computed stress intensity factors can be compared. 
However, the computed values of K x can be compared to those obtained from the finite element 
analysis of Miller, et al. (1992) of a similar panel test. The comparison is described in Potyondy (1993) 
and Potyondy, et al. (1994b) where it is shown that the trends and values from both analyses are similar. 


Fatigue Life Prediction. A life prediction has been made using the computed stress intensity 
factor history and a Paris fatigue crack growth model. Fatigue crack growth data from a flat 2024-T3 
panel tested at stress intensity factors in the range of 50-80 ksi-sqrt (in) was made available to the 
authors by Boeing [Worden (1993)]. The crack growth equation is of the form 

(10 ’ 000) || = ( ; % !1 J <3) 
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and for 2024-T3 aluminum T-L crack growth ( R = 0 and ambient lab environment) 

C = 28 

n = 4 9 * forjK max (ksiVinch) ( 4 ) 

By using the computed values of and assuming that K, varies linearly between the data points in 
igure 10, life predictions were made and are compared with measured life in Figure 1 1 (K „ is defined 
below by Eq. (6); for a = 0,K^= K,). 


The total predicted life is obtained by integrating the crack growth model until K, reaches K 
the apparent toughness of the material. K app is a function of thickness and crack orientation relative to 
material rolling direction. The values of K app for this material and thickness are 


K-app — 100 ksiVinch (T-L crack growth) 

K-app ~ 105 — 120 ksiVinch (L - T crack growth) 


C u Ck a i igned witl l the roUin § direction, while L-T crack growth refers to 
P®^. endl ? ular to the rolhng direction. For the panel tested, the rolling direction is the 
g u ina direction; thus, cracks growing from the initial saw cut are experiencing T-L crack growth. 

consfo^hl ted Tf- 116 ° f Ki 1° T L ^ B = °' 77 iS 112 ksi ' s( l rt 0 n ). At this point, the crack has curved 
th l erab y ^ r 1 l onented at nearl y 45 degrees with respect to the two orthogonal rolling directions; 
thus, the value of K app will lie within the range of T-L and L-T crack growth. For the purposes of 

making a total life prediction, it is assumed that K, = K app for L/B = 0.77. 


well wi T he Ufe Potion for the vahdation problem compares 
I th mea sured life. For all crack lengths considered in the simulation, the predicted life is never 
Tn°th e than w P ercent greater than the measured life. In addition, the total life prediction is 628 cycles 

SS s r gr0Wth ceas . ed after 626 c y cles because of the large flap^hat haJ fSmS in tire 
with further pressurization impossible. Thus, the total predicted life is in excellent agreement 

with the measured life; however, such an extremely close comparison is surely serendipitous 


Life predictions within a factor of two are considered good for this test for the following reasons: 

The crack growth rates throughout the test are high; thus, the predicted life is sensitive 
to small changes in the computed values of because of the large coefficient of 4 9 
in the Paris crack growth model 


During the final stages of crack growth, the computed values of Kj are approaching 
ap P ’ thus, the crack growth is about to become unstable. It is likely that large plastic 

iS Were Present at the crack tips, thereby casting doubt on applicability of the 
LE ^ assumptions. However, this effect is offset by the use of Paris model 
coefficients that were obtained for stress intensity factors above 50 ksi-sqrt (in) 


A final observation is that the prediction of slower crack growth than observed supports the 
ypothesis that the presence of k z will increase the crack growth rates from those obtained by 
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considering only K I . It is expected that the experimental work of Zehnder, et al. (1992) and Viz and 
Zehnder (1994) will shed more light on this issue and may result in a modified procedure for 
determining life. In lieu of such a procedure, a simple empirical procedure for incorporating and 
bounding the effect of the additional stress intensity factors on the crack growth rate is proposed here. It 
assumes that a fraction of the energy from the additional stress intensity factors contributes to fatigue 
crack growth. Replace the value of Kj in the fatigue crack growth model with an effective stress 
intensity factor K^, that is related to an effective energy release rate G^, via 


G e/r =G 1 + a(G 2 +G 3 +G 4 +G 5 ) ; 0<a<1.0 

where E is Young’s modulus, G t are the separate components of the energy release rate that are related 
to the four stress intensity factors by Eq. (1), and a is an empirical coefficient that scales the 
contribution from the additional stress intensity factors. 


The empirical procedure is applied to the validation problem, and the resulting life predictions 
for a equal to 0.0 and 1.0 (expected lower and upper bounds) are compared with measured life in Figure 
1 1. The a coefficient shifts each computed data point to the left by a fixed amount, e.g., the last 
computed data point will coincide with the measured value by setting oc = 0.4 . Also, total life is 
decreased as a is increased. During crack growth, although the ratio of k^/Kj increases by almost 40 
percent, the contribution of K I to the total energy release rate decreases by less than 4 percent 
suggesting that may have a minor effect on crack growth for this class of problems. 


Discussion 


The present analysis methodology allows the cracks to curve and accounts for the nonlinear 
bulge-out effect and the presence of stiffening elements by using a geometrically nonlinear solution 
procedure. Also, the interaction between the mechanical loads acting on the superstructure and the local 
structural response near the crack tips is accounted for by employing a hierarchical modeling strategy. 
The computed stress intensity factors exhibit the expected trends and their values are reasonable. The 
computed crack trajectory and predicted life agree well with the measured quantities from a full-scale 
pressurized panel test. The effect of the additional stress intensity factors, in addition to Kj , on crack 
growth rate can be incorporated into a life prediction using an effective stress intensity factor. 


The unsymmetric, out-of-plane deformation experienced by a crack in a pressurized fuselage 
structure can be seen clearly in the geometrically nonlinear finite element models (see Figures 5 and 8) 
and can be accounted for by the computation of . However, the effect of k% on crack growth rate is 
still uncertain. It is expected that further testing by Zehnder, et al. (1992) and Viz and Zehnder (1994) 
will determine this effect. The slight underprediction of crack growth rate for the validation problem 
may indicate that the presence of increases the crack growth rate beyond that obtained by only 
considering K : ; however, this result is preliminary. Other factors, such as the lack of applying a 
longitudinal loading to the global shell model, may also account for the slight underprediction of crack 
growth rates. Further comparisons between predicted and measured lives must be made before more 
definitive conclusions can be drawn. 
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tr ■ ,J he ™™ tangential stress crack growth directional criterion used here predicts a crack 
m ^good agreement with the measured trajectory from a full-scale pressurized panel test (see 
f ^ impending upon the accuracy required in the trajectory determination, this simple criterion 

may be sufficient for most practical situations. v ° 


The fact that the boundary conditions on the 2x2 bay stiffened panel model were not altered 
during crack growth; i.e„ the entek was only grown in the 2x2 bay model, and notTeriousWfect 
mfnbSh 3 ? ° f th ® 1 trajectory and life predictions. The kinematic boundary conditions are affected 

mally by crack growth near the center of the panel and must only be updated when the crack begins 
l a P pr °T h rS he pane J bpuiKtary. The crack behavior is affected to a much greater extent by the live g 
p essure loading, and the effect of this loading is being recomputed after each step of crack growth. 


Restarting the geometrically nonlinear finite element computation from zero load after each 
^ a significant computational effort. The computational efficiency could 

be improved by using the previously computed equilibrium state as an initial condition for the 

anH^nDr^^nono? 10 !? after . the crat * has grown. The procedure is described in the paper of Potyondy 
growtlfsimu?ation sign^ficamly 016 "* 3 ^ 011 rcdUCe the ^ com P ut£ " io ” P^orm a crack 


that fastene F flexibility alters the stress distribution locally in the fuselage skin which 
influences the behavior of cracks in these regions. The effect of neglecting fastener flexibility in the 

sh^flH t S , H y t S1S me * oc J olog y f^ould be examined. Also, the repeatability of the panel test experiments 

oald ^ d — ed along with the effect of small variabilities, such as fit-up stresses and rivet forces, 
on crack growth behavior. 


CONCLUSIONS 


The present methodology provides an effective engineering approach to simulate crack growth in 

^ here crack ™ not know a priori. It accounts for the dcSrinant 

SjS 5 , aff ^ ctl . n 8 ^rack growth behavior (the bulge-out effect; presence of stiffening elements; 
fhmi^dth ‘ a 2? i™? add ; tlon to the internal pressure, acting on the superstructure) and integrates 
craS rr M ^^P^eters to support trajectory and life predictions. The predicted 

J 2 d bf a S re ® well with the measured quantities from a full-scale pressurized panel test 
dicatmg the effectiveness of the methodology. A logical extension of the present methodology would 
be to account for plasticity effects at the crack tips, as these effects may become significant fo/long 
cracks encountered during residual strength computations. 8 g 
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fracture simulation program analysis program 


Figure 1: Discrete crack growth analysis methodology. 




bending loading, 



Anti-symmetric 
bending loading, 


The four loading modes and their corresponding stress intensity factors 
for a through crack in a thin plate idealized using plane stress elasticity 
to account for the membrane loads and Kirchhoff plate theory to account 
for the bending loads. 


Figure 2: 











2x2 bay stiffened panel model 


Figure 3: The three hierarchical modeling levels. Each model obtains its 

hierarchy boUndary condltlons from the proceeding model in the 
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tear straps 



Figure 4: The 2x2 bay stiffened panel model as represented within FRANC3D as a 

geometrical model. The structural components and their cross sections 
are labelled. 
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Figure 5: Deformed structure during crack growth for total crack lengths of 6, 8, 

10, and 12 inches for the trial problem. Magnification factor for all 
images is five. 
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Figure 6: Computed crack growth trajectory for the trial problem. The black 

circles denote computed crack tip locations. 



Relative crack length, LjB 


Figure 7: Computed stress intensity factors at right crack tip for the trial problem. 
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Figure 8: Deformed structure during crack growth for total crack lengths of 5, 7, 9 

and 1 1 inches for the validation problem. Magnification factor for all 
images is five. 
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Stress Intensity Factor 


Frame Midbay Frame 



Figure 9: Comparison between computed and measured crack growth trajectories 

for the validation problem. The black circles and open boxes denote 

computed and measured crack tip locations, respectively. 



Relative crack length, L/B 


Figure 10: Computed stress intensity factors at right crack tip for the validation 

problem. 
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Figure 11: Comparison between predicted and measured fatigue life. The life 

prediction is made using the effective stress intensity factor of Eq. (6). 




